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We study beams of medium-energy electrons and positrons which obey the Dirac equation and scatter from nuclei with 

Z > 100. At small distances the potential is modelled to be that of a charged sphere. A large peak is found in the probability 

of positron penetration to the origin for Z ~ 184. This may be understood as an example of Klein tunnelling through the 

Coulomb barrier: it is the analogue of the Klein Paradox for the Coulomb potential. 
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I. QUANTUM TUNNELLING AND KLEIN 
TUNNELLING 

One of the principal characteristics of quantum me- 
chanics as opposed to classical mechanics is that particles 
of energy E can pass through regions of a potential V(r) 
in which their kinetic energy T = E— V is negative, albeit 
with amplitudes which are exponentially suppressed. In 
two pioneering papers 70 years ago which demonstrated 
the importance of quantum tunnelling in the understand- 
ing of nuclear decays, the probability of an a-particle to 
tunnel from the potential well of the nucleus through the 
Coulomb barrier was calculated 0. The calculation in 
the quasi-classical case is straightforward || : the transi- 
tion probability for decay w is given by 



ir - cxp \ -2 J ^J2m [V(r) - E]j = cxp {-2irZam/p} 

(1) 

for V(r) = Za/r where r c = Za/E is the classical turn- 
ing point, and m and p are the mass and momentum of 
the a-particle. The integral in the exponent is over the 
classically forbidden region in which T < 0. Similarly a 
positron of energy E, momentum p, and mass m incident 
on a nucleus of charge Z can reach distances smaller than 
r c . If p = |^(0)| / \ip(0)\ e i is the ratio of the probabil- 
ity of a positron penetrating a Coulomb barrier to reach 
the origin compared with the probability of an electron 
of the same energy, then essentially the same calculation 
gives non-relativistically 



p = e- 2ny 



(2) 



where y = Zam/p so that p decreases exponentially with 
Z for fixed p. If the particles are relativistic and satisfy 
the Dirac equation then Eq. is still obtained provided 
now y = ZaE/p [El. For relativistic problems the ki- 
netic energy T = E — V — m = y m 2 + p 2 — V — m 



so that, classically, T is still positive definite but now 
r c = Za/(E -m). 

At the same time as the papers showing the effect of 
quantum tunnelling in a-decay, Klein Q showed that 
electrons in the Dirac equation could in principle tunnel 
through a high repulsive barrier. This is the famous Klein 
Paradox. Calogeracos and one of us (ND) have recently 
reassessed this phenomenon |6) and call such quantum 
tunnelling without the expected exponential suppression 
Klein tunnelling. To obtain Klein tunnelling it is essential 
that hole states (corresponding to negative energy states 
for free particles) as well as particle states are considered 
and allowed to propagate. For a particle of momentum p 
and energy E incident on a Klein step of height V, a hole 
state of momentum — q will propagate under the barrier 
provided 



V = \J m 2 + p 2 + \J m 2 + q 2 , 

which is possible for V > 2m. In terms of the particle 
kinetic energy under the barrier 



T = E-V 



— y m 2 + q 2 < —2m, 



and where T < —2m hole states can propagate without 
exponential suppression. T < —2m for a Coulomb po- 
tential corresponds to penetrating under the barrier to 
distances r < tk, where tk = Za/(E + m) < r c is the 
Klein distance. While tunnelling from r c to tk is expo- 
nentially suppressed, tunnelling from tk to r — is not; 
indeed the amplitude may even be enhanced. This is be- 
cause the effective potential in a relativistic theory with a 
potential V which is the time component of a four- vector 
is given by 



2mV eff (r) = 2EV(r) - V 2 (r). 



(3) 



In the region < r < tk, V( r ) > E + m > E and 
so the effective Coulomb force is attractive || when r 
is small enough [note that the force becomes attractive 
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when dV e ff/dr changes sign, not when V e ff (r) changes 
sign]. 

When the Coulomb potential is investigated in the 
Dirac equation, the ground state energy E\ = at 
Za = 1 and becomes complex for Za > 1 1, This is 
a consequence of the V 2 {r) term in Eq. (||) which leads 
to the "collapse" of the particle to the origin, and hence 
to a problem which is not well-defined . So the theory 
breaks down at .Z max = 1/ot — 137. Since superheavy 
nuclei can be constructed in heavy ion collisions with 
values of Z larger than 137, this limitation on Z cannot 
be physical. A modified Coulomb potential which takes 
account of the finite size of the nucleus must therefore 
be used so that the singularity at small r is smoothed 
out. When this is done there seems to be no restriction 
on Z; furthermore bound state energies become negative 
for sufficiently large Z. These arc interpreted as bound 
positron states ||. For Klein tunnelling through the 
(modified) Coulomb potential it was conjectured Q that 
it might occur at values of Z large enough to obtain nega- 
tive energy bound states. For the simplest such potential, 
due to Pieper and Greiner 63], the first negative energy 
bound state occurs at Z — 147. Hence the prediction of 
Calogeracos and Dombey || that Klein tunnelling may 
occur for nuclei with Z m 150 or above. If it occurs, Klein 
tunnelling would lead to a breakdown of the exponential 
suppression given by Eq. (|J) . 

In the Dirac equation in one-dimension, transmission 
resonances (with zero reflection coefficient) |l(| have been 
demonstrated for electron scattering off square barriers 
[11 1 and off smooth potential barriers of Woods-Saxon 
[12 1 or Gaussian jl3) form. These provide examples of 
Klein tunnelling. In each case, the transmission reso- 
nances occur when the corresponding attractive potential 
becomes supercritical (see next paragraph). 

Since it is unlikely that superheavy nuclei with Z > 
150 can be prepared and stay around for long enough 
for positrons to be scattered off them, we study positron 
scattering off nuclei with Z > 150 numerically by solving 
the Dirac equation for a modified Coulomb potential. We 
thus attempt to repeat Klein's analysis for a modified 
Coulomb potential in place of a potential step. We do 
indeed find that positrons scattering off nuclei with Z > 
150 and especially with Z > 170 no longer satisfy Eq.(^|). 
We find, in addition, an extremely sharp peak in p near 
Z = 184. We now outline the calculation. 



II. MODIFIED COULOMB POTENTIAL: BOUND 
AND CONTINUUM STATES 

Following Pieper and Greiner ||, the modified 
Coulomb potential experienced by an electron or positron 
beam is taken to be the potential arising from a homoge- 
neously charged sphere of radius b. More explicitly, the 
time component V(r) of the 4-vector potential is given 
by 



V(r) = ±^f(r/b), 

where the dimensionless potential shape function / is de- 
fined as 

/(r ./ 6 ) = -| + l(r) 2 r <b 



f{r/b) 



r > b 



and b — (1.2)^4.3 fermi. We let ipi and tp2 be the 'large 
and small' radial wave functions used to construct the 
Dirac spinor corresponding to a total angular momentum 
of j. We employ the variables r = ±1, and k = j + 
i, so that the parity P of the spinor is given by P = 
(— iy + i — ±1. In the notation of Rose Q], with k = rk, 
the coupled radial equations may be written 



^2 



rk 



'-ip 2 = (m + V - E)i/)i 



ip'i H Vi = (m-V + E)ip 2 - 

r 



(4) 



(5) 



There are two distinct problems to be considered: (a) 
bound states; and (b) continuum states. For the bound 
states, we adopt the boundary conditions and normaliza- 
tion given by: 



iMo) = V2(o) = o, 



(^(r) + V 2 2 (r))rfr = l. 



We now choose j and t, integrate out from the origin, 
and search for those energies which lead to spinors that 
have the desired number n\ of nodes in the 'large' ra- 
dial component tpi, and vanish at large distances. If 
the quantity t = j + represents the orbital angular- 
momentum quantum number in the first two compo- 
nents of the Dirac spinor, and v = (ni + 1 + t), then 
the usual spectroscopic designation is written v£j , where 
I = {0, 1,2,...}- {s, p, d,.. .}. The j = |, I = states, 
for example, correspond to t = —1. We calculate the first 
few energy eigenvalues as a function of Z and our results 
agree with those given by Pieper and Greiner in their Ta- 
ble 3. In particular, we have found the critical values of 
the atomic number Z that yield the energies E = (the 
threshold for bound positrons) and the supercritical val- 
ues of Z corresponding to energy E = —m (correspond- 
ing to spontaneous positron production) ; some results are 
presented in Table 1 





E = 


E = —m 


lSl/2 


146.7 


170.4 


2«l/2 


195.7 


237.0 


2^1/2 


168.1 


183.8 


3S1/2 


260.1 


316.5 



Table 1: Values of Z for which E — and E = —m 
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In the case of the continuum states, we adopt the follow- 
ing method. At very small distances r << 6, we assume 
that the large and small radial functions have the asymp- 
totic pure-power forms CV 3 and the ratios given (after 
some elementary corrections) by Rose 0. In this region 
there is then only one free parameter, which we take to 
be the amplitude C\ e , or, for positrons, Ci p , of the large 
component ipi (r); the small component is smaller here by 
the factor r 2 . We fix j, t, and the incident momentum p, 
and we integrate outwards, to a point well beyond the 
classical turning point r c . At these distances the radial 
components have the large-r Coulombic asymptotic form, 
which is approximately sinusoidal with ^i-amplitude, say 
A%; in this region, the components are also (almost) ex- 
actly out of phase, so that when ipi = Ai, we have tp2 = 0. 
We now search, in each case, for the value of C\ so that 
this asymptotic amplitude has the value A\ = 1. Since 
the numerical process is, in principle, direction invariant, 
this corresponds to an 'experiment' in which the final 
amplitude of an incoming beam with unit amplitude is 
determined near the origin. Since we are calculating wave 
functions at the origin, only I = states contribute and 
thus we restrict the analysis to r = —1 and k = 1. For a 
given Z we can express the positron-electron ratio in the 
form p = [Cip/Cie] 2 . 



III. THE PEAK 

As we are not interested in the detailed behaviour of 
the nuclear charge distribution we require b << tk- We 
also require r c » r# so that the normal exponential 
suppression will be obtained for Z not too large. We en- 
sure that this is so by choosing the momentum p in the 
range 0.1m — 0.4m. For this momentum range and with 
Z < 150 we expect normal exponential suppression of the 
positrons by the modified Coulomb potential according 
to Eq.(j^): this is demonstrated in Fig. 1 where our calcu- 
lation of p for Z — 100 is close to e~ 27ry asymptotically, 
showing that our method of calculation is satisfactory. 

We now look at C\ p /C\ e for values of Z > 140. The 
principal discovery we report here is the existence of a 
pronounced peak in p — [Ci p /Ci e ] 2 very near Z = 184. 
This peak is exhibited in Fig. 2 for the case p = (0.4)m. 
At the peak maximum p rj 2. The peak is so sharp 
that it seems to correspond to a singularity of some sort. 
It is convenient to look at the deviation from exponen- 
tial suppression by introducing a Klein " logarithmic form 
factor" R so that \C\ v jC\ e \ = e R e~ 7ry , where R is given 
explicitly by 

R(z, P ) = in \c lp /c le \ + *y = \ Hp) + *v- (6) 

Thus R constant corresponds to normal exponential sup- 
pression. In Fig. 3 we exhibit more detail near the peak 
by plotting R in terms of Z for the two cases p — 0.1m 
and 0.4m. The peak shifts slightly with momentum, with 
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FIG. 1. A comparison of p (solid) with e 2,ry (dashed) for 
Z = 100. 



ln(p) 




I i i i i I 



140 150 160 170 180 190 



FIG. 2. ln(p) as a function of Z 



a limit, for p — > 0, very close to the p = 0.1 position 
shown. More specifically, we find peaks in p(Z,p) at posi- 
tions Z p given by Z 0A = 184.8 and Zq.i = Z .o2 = 183.8. 
This is fully consistent with a peak in p(Z,p = 0) at Zq = 
183.8. 



IV. DISCUSSION 

The enormous peak near Z = 184 in the scattering 
of positrons by superheavy nuclei was unexpected: it 
involves an enhancement in p by a factor 10 8 between 
Z = 179 and Z = 184. While a demonstration of a de- 
viation from the pure exponential suppression of Eq.(|^) 
resulting from Klein tunnelling was the purpose of this 
calculation, we initially expected it to be a simple max- 
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FIG. 3. R as a function of Z for p — 0.4m (dashed) and 
p = 0.1m (solid) 

imum in p at the first supercritical value of Z = 170.4 
as is the case with the one-dimensional examples referred 
to above, where supercriticality and the maximum of the 
transmission coefficient coincide. While the peak that we 
have found does not correspond to the first lsi/2 super- 
critical state, it does correspond to the second 2pi/ 2 su- 
percritical state (see Table 1). We show in a forthcoming 
paper [fL4| that this correspondence exists independently 
of the exact form of the potential. This is also in agree- 
ment with the result of Muller, Rafelski and Greiner ]l5|| 
who predicted a resonance in the negative energy contin- 
uum in the S1/2 state at Z = 184. 

To explain why the peak occurs in the state which cor- 
responds to the 2pi/2 bound state of the electron is an 
exercise in Dirac's hole theory. The vacuum state for 
our modified Coulomb potential at supercriticality con- 
tains a vacant 2p 1 / 2 bound state of an electron of energy 
E = — m; this in hole theory represents a positron state of 
zero kinetic energy. But to compare the original electron 
state with a positron state requires charge conjugation: 
equations (U^) are invariant under the combined oper- 
ation E — > -E;V — * -V\t — > —T;ipi 45 tp2- So the 
vacant electron 2p 1 / 2 bound state corresponds to a reso- 
nant S\/2 positron scattering state since j = | for both 
states but r changes sign.Q 

The (vacant) electron 2pi/ 2 bound state of zero ki- 
netic energy at supercriticality can thus be considered 
as a positron resonance of zero kinetic energy. In one 
dimension non-relativistically a zero energy resonance is 
a transmission resonance where the transmission coeffi- 
cient is unity and the reflection coefficient is zero jl6| . It 



1 We thank Dr. X. Artru for pointing this out. 



seems likely that we have found a similar effect here in 
three dimensions since our peak occurs as close to zero 
kinetic energy (corresponding to the supercritical energy) 
as we can calculate. We should emphasise that we calcu- 
late the quantity p which is not a scattering cross section. 
Neither is it strictly a transmission coefficient. Neverthe- 
less our peak is more like a transmission resonance than 
a resonant cross section. We will examine this point in 
more detail in a forthcoming paper. 

We note finally that our example here of positron scat- 
tering by a modified Coulomb potential shows that the 
phenomenon discovered by Klein of Dirac particles tun- 
nelling through repulsive potentials is a general charac- 
teristic of the Dirac equation in the presence of strong 
fields. 
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